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THESIS  ABSTRACT 

SOLID-PROPELLANT  GENERAL  DEFORMATIONS  AND 
THEIR  EFFECTS  ON  MOTOR  PERFORMANCE 

David  Franklin  Smith 

Master  of  Science,  December  9»  19?6 
(B.A.E. , Auburn  University,  1975/ 

48  Typed  Pages 

Directed  by  Richard  H.  Sforzini 

A simplified  analysis  of  the  propellant  strains  in  solid-pro- 
pellant rocket  motors  with  circular  perforated  grains  is  presented. 

The  analysis  extends  the  work  of  Vandenkerckhove  to  include  flexible 
motor  cases.  Each  of  the  propellant  strains  is  assumed  independent 
of  axial  position  along  the  propellant  grain.  Strains  calculated  from 
this  analysis  are  compared  with  the  axially  varying  strains  obtained 
from  the  finite  element  analysis  of  Brisbane.  The  axial  and  radial 
strain  comparisons  show  poor  agreement  for  the  strains  predicted  by 
this  theory,  however,  the  tangential  strain  correlation  is  excellent 
over  a large  portion  of  the  propellent  grain. 

This  analysis  for  tangential  propellant  strain  is  used  to  de- 
termine general  propellant  deformations  and  their  effects  on  rocket 
motor  internal  ballistics  performance  predictions.  Performance  pre- 
dictions with  these  deformation  effects  are  compared  with  experimental 
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data  on  two  rocket  motors.  In  both  cases,  the  performance  predictions 
with  the  propellant  deformation  effects  are  closer  to  actual,  motor 
performance  than  are  the  predictions  without  the  deformation  effects. 
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I.  INTRODUCTION 


One  of  the  problems  in  the  prediction  of  solid  propellant  rocket 
motor  performance  is  the  variation  between  the  actual  rocket  motor  and 
laboratory  propellant  burning  rates.  In  the  laboratory,  a sample  of 
the  solid  propellant  to  be  used  in  a given  full  scale  rocket  motor  is 
burned  to  determine  the  "exact"  burning  rate.  This  is  accomplished  by 
firing  a small  ballistic  test  motor  loaded  with  propellant  or  by  burn- 
ing a thin  strand  of  propellant  under  laboratory  conditions.  Neither 
the  ballistic  test  motor  nor  the  strand  burner  technique  gives  consist- 
ently reliable  results.  The  apparent  burning  rate  of  the  propellant  in 
the  actual  motor  may  be  higher  by  1C$  or  more  than  that  indicated  by  the 
strand  burner  or  ballistic  test  motor. 

One  explanation  of  a portion  of  this  burning  rate  discrepancy  lies 
in  the  deformation  of  the  solid  propellant  during  motor  operation.  The 
circular  perforated  portion  of  the  propellant  grain  in  a solid  propellant 
rocket  motor  (SRM)  will  deform  due  to  the  pressurization  from  the  com- 
bustion chamber  gases  when  the  motor  is  fired.  This  deformation  will  re- 
sult in  a larger  burning  perimeter,  and  thus  in  a higher  apparent  burn- 
ing rate.  The  actual  burning  rate  of  the  propellant  referred  to  a given 
pressure  may  not  vary  between  the  strand  burner  or  ballistic  test  motor 
and  the  actual  SRM,  but  the  increase  in  burning  perimeter  will  result 
in  more  propellant  being  burned  in  a given  time  interval  than  expected. 

Thus,  over  the  total  burning  time  of  the  propellant,  the  bum.ng  rate 
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seems  to  be  higher  than  predicted.  This  variation  in  the  apparent  bum- 
I ing  rate  is  directly  proportional  to  the  amount  of  propellant  defor- 

■ mation  at  the  bore  of  the  grain. 

The  burning  rate  discrepancy  is  especially  important  on  new  motor 
I designs  where  the  "scale  factor"  on  burning  rate  is  unknown  or  un- 

certain. Until  one  or  more  full  scale  motors  are  fired,  actual  motor 
I performance  is  seldom  accurately  predicted  based  on  laboratory  burning 

| rate  information.  Thus,  it  is  important  to  isolate  and  determine  the 

factors  that  affect  the  burning  rate  determination. 

I In  an  SRM  of  the  Space  Shuttle  type  with  a propellant  web  thick- 

ness of  40  inches,  the  initial  burning  bore  perimeter  may  increase  by  as 
I much  as  k%  due  to  propellant  deformation,  with  a corresponding  increase 

■ in  the  apparent  burning  rate.  This  example  is  not  a limiting  case  and 
other  rocket  motors  may  show  even  greater  deformations. 

I Most  SRM  performance  prediction  techniques  utilize  a software  com- 

puter program  in  the  design  stages.  The  ability  to  incorporate  the  pro- 
I pellant  deformation  effects  in  the  design  of  a SRM  can  improve  the  per- 

■ formance  predictions  on  the  new  motor  and  perhaps  point  out  some  of  the 
prediction  errors  on  older  motors. 

I It  is  the  objective  of  this  Investigation  to  present  a simplified 

analysis  of  the  propellant  deformation  in  SRM's.  This  investigation  is 
1 confined  to  the  deformations  of  circular  perforated  (c.p.)  grains  since 

^ many  propellant  grains  are  all  or  mostly  c.p.  In  addition,  due  to  the 

simplicity  of  the  geometry,  the  c.p.  grain  is  a good  first  test  for 
J the  hypothesis  with  possible  future  efforts  aimed  at  a strain  analysis 

of  the  more  complicated  star  or  wagon  wheel  grain  configurations.  : 

I ! 
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The  grain  deformations  predicted  by  this  analysis  are  due  to 
uniform  pressure  loading.  Also,  this  analysis  predicts  the  general  def- 
ormations of  the  grain,  as  opposed  to  local  deformations  along  the  axis 
of  the  grain.  Local  deformations  can  be  caused  by  end  effects,  pressure 
gradients  in  the  combustion  chamber  or  flight  acceleration1’2  and  are 
not  considered  in  this  analysis. 

In  the  present  analysis , the  method  of  Vandenkerckhove  2 for  rigid 
motor  cases  has  been  modified  to  include  the  effects  of  flexible  cases. 
In  addition,  an  existing  SRM  performance  prediction  program  is  utilized 
to  show  the  effect  of  these  propellant  deformations  on  the  performance 
predictions  for  two  motors.  Experimental  data  on  these  two  motors  is 
compared  with  the  performance  predictions  with  and  without  the  propel- 
lant deformation  effects.  The  increased  accuracy  of  the  prediction  with 
the  deformation  effects  is  evidence  of  the  applicability  of  this  inves- 
tigation. 


II.  ANALYSIS 


A tubular  grain  solid  propellant  rocket  motor  consists  of  a 
cylindrical  body  of  propellant  which  has  a circular  perforation  along 
the  centerline  (See  Figure  1)  and  is  enclosed  by  a cylindrical  case. 

This  type  of  solid  propellant  grain  is  usually  referred  to  as  a circu- 
lar perforated  or  c.p.  grain.  When  analyzing  the  stresses  and  strains 
in  such  a motor,  it  is  convenient  to  model  the  motor  as  a body  of  pro- 
pellant enclosed  by  a thin  walled  flexible  case. 

A section  of  the  c.p.  grain  is  shown  in  Figure  2 with  the  deformed 
and  undeformed  propellant  surfaces  as  indicated.  The  deformation  is 
greatly  exaggerated  for  clarity,  but  the  burning  perimeter  (s)  of  the 
propellant  is  seen  to  Increase  by  S€g  where  6^  is  the  tangential  strain 
in  the  propellant  at  the  bore  of  the  grain.  The  propellant  deformation 
is  symmetrical  with  respect  to  the  axis  along  the  bore  of  the  c.p.  grain 
and  consists  of  a radial  displacement  of  the  grain  that  is  constant 
along  the  circumference  of  the  cylinder  but  varies  along  the  radius. 

That  is,  the  deformation  is  a function  of  radial  position  in  the  grain 

3 

but  at  a given  radius  the  tangential  strain  is  a constant  . Also,  each 
of  the  propellant  strains  is  assumed  to  be  indeper^ent  of  the  axial  po- 
sition. The  validity  of  this  assumption  will  be  tested  when  the  com- 

4 

parlson  with  the  finite  element  theory  of  Brisbane  is  made. 
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Figure  Z.  Cross  Section  Segment  of  a Solid  Propellant 
Circular  Perforated  Grain  with  and  without 
Pressurization 
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One  may  assume  that  the  radial  stress  (tf  ) is  essentially  zero  in 
the  thin  walled  motor  case  and  from  the  hoop  stress  formula  the  tangen- 
tial stress  in  the  motor  case  is  given  by 


°e  - VA  « 

where  Rg  is  the  radius  of  curvature  of  the  motor  case  (i.e.  the  outside 
radius  of  the  c.p.  grain),  X is  the  case  thickness  and  P is  the  radial 
stress  in  the  propellant  at  the  case  wall. 

For  a non-rigid  motor  case,  the  tangential,  radial  and  axial  strains 
are  non- zero  and  are  given  by  Hooke's  Law  as 

€e  = - vfar  + tfz)]/E  * (2a) 


er  *=  Or  - v(tf0  + tfz)]/E  ¥ 0, 
€z  " Oz  " v(*0  + tfr)]/E  * 0. 


(2b) 


(2c) 


In  Equations  (2a)  through  (2c),  the  subscripts  0,  r and  z refer  to 
tangential,  radial  and  axial  directions,  respectively.  The  strains  are 
denoted  by  €,  the  stresses  by  cT , the  elastic  modulus  by  E and  Poisson's 
ratio  by  v. 

If  the  rocket  motor  were  a closed  pressure  vessel,  then  the  axial 
stress  in  the  case  would  simply  be  one  half  of  the  hoop  stress  (See  Fig- 
ure 3)1  however,  in  a rocket  motor  the  situation  is  not  as  simple. 

There  is  a uniform  pressure  force  on  the  head  end  of  the  motor  and  a 
non-uniform  pressure  force  on  the  case  shoulders  forward  of  the  throat 
and  in  the  exit  cone  of  the  nozzle.  For  the  purpose  of  a static  anal- 
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ysis,  there  are  also  reverse  effective  forces  acting  through  the  center 
of  gravity  of  the  body  equal  in  magnitude  to  the  motor  thrust  (t)  times 
the  mass  fraction  (m^)  of  that  portion  of  the  body.  In  Figure  4,  M in- 
dicates mass  and  m^,  indicates  the  mass  fraction  which  is  the  mass  of  a 
portion  of  the  body  divided  by  the  total  mass  of  the  body.  The  sub- 
scripts i,  2,  3 and  s indicate  the  payload,  forward  portion  of  the  body, 
aft  portion  of  the  body  and  the  entire  vehicle,  respectively.  For  exam- 
ple, the  reverse  effective  force  acting  on  the  thrust  skirt  is  equal  in 
magnitude  to  the  mass  fraction  of  the  payload  (m„  ) times  the  motor 
thrust.  It  Is  obvious  that  a rigorous  analysis  of  the  axial  stress  in 
the  motor  case  becomes  very  complicated,  since  this  stress  is  a function 
of  both  time  and  position  along  the  axis  of  the  motor  (i.e.  the  mass 
fraction  of  the  free  body  will  be  different  at  every  axial  location, 
and  the  pressure  and  thrust  may  both  be  functions  of  time). 

It  is  desirable  to  choose  a constant  axial  case  stress  for  the  pur- 
pose of  simplifying  this  analysis  while  observing  that  the  case  strain 
than  most  greatly  affects  the  propellant  bore  strain  will  be  the  tan- 
gential motor  case  strain.  The  propellant  has  essentially  a constant 
volume,  since  Poisson's  ratio  is  close  to  0,5  for  most  composite  solid 
propellants.  Therefore,  a tangential  case  strain  that  causes  a small 
change  in  case  diameter,  and  thus  a small  change  in  the  outside  diam- 
eter of  the  case-bonded  grain,  will  produce  a somewhat  larger  increase 
in  the  grain  bore  diameter.  The  axial  case  strain  not  only  has  a small- 
er effect  on  the  grain  bore  deformation  for  a given  strain,  it  is  also 
smaller  in  magnitude  than  the  tangential  case  strain.  Thus,  the  valid- 
ity of  any  assumption  on  the  axial  case  stress  must  be  considered  with 
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respect  to  the  effect  on  the  tangential  motor  case  strain.  It  is  also 
important  to  note  that  the  effect  on  the  propellant  deformation  at  the 
bore  of  the  grain  from  any  of  the  motor  case  strains  will  be  small  com- 
pared to  the  deformations  due  to  the  internal  pressure  loading. 

The  minimum  axial  case  stress  (assuming  that  the  motor  case  is  al- 
ways in  tension)  would  be  zero  and  a clear  maximum  would  be  one  half  of 
the  hoop  stress  (since  this  corresponds  to  a closed  pressure  vessel). 

An  application  of  Hooke's  Law  for  the  situation  where  the  axial  motor 
case  stress  is  zero  and  with  a nominal  value  of  0.3  for  v of  the  case 
yields  (for  d ^ •=  d^  = 0) 


and 


E € 
c zc 


-0.3  d 


0c 


E 

c 0c 


0c' 


(3a) 

(3b) 


where  the  c subscript  refers  to  the  motor  case.  For  a maximum  axial 
case  stress  of  one  half  of  the  hoop  stress,  we  have  (for  d ■=  0.5  da 

ZC  w 

and  d ■=  0) 
rc 


Vsc  ■ °-5  S6o  - °-3(c,ec)  ■ °-2  ff6c  (3c) 

and 

eoe6c  • aec  - tfec>  ■ °-85  °eo-  <*> 

Comparing  Equations  (3b)  and  (3d),  one  sees  that  the  tangential 
case  strain  changes  only  15#  between  the  maximum  and  minimum  assumed 
values  for  the  axial  case  stress.  We  thus  have  a small  change  in  a 
quantity  that  has  a small  effect  on  the  propellant  deformation.  In  ad- 
dition, the  maximum  axial  case  stress  assumed  in  this  argument  is  clearly 


] 
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greater  than  any  real  axial  case  stress.  Thus,  by  assuming  a zero  axial 
motor  case  stress,  one  nay  expect  an  even  smaller  change  in  the  tangen- 
tial case  strain  and  a resulting  minor  difference  in  the  actual  propel- 
lant deformation  at  the  bore  of  the  grain.  It  should  also  be  noted,  as 
shown  in  Figure  4,  that  for  hi$i  acceleration  vehicles  (i.e.  hif£i  thrust 
to  weight  ratio)  the  axial  case  stress  may  well  be  close  to  zero. 

For  the  thin  walled  case,  with  the  radial  and  axial  stresses  assum- 
ed to  be  zero,  Equations  (2a)  through  (2c)  are  rewritten  as 

feo  ■ W 

- ‘VeA-  ('•*>) 

and 

«»  ■ -VeA-  <*>> 

With  Equation  (l),  one  may  further  modify  Equations  (4)  to  become 

€0c  ■ <*> 

eic--VoB.V<Vc)'  <5b> 

and 

Assuming  that  the  propellant  is  bonded  to  the  case,  the  strains  at 
the  outer  radius  of  the  propellant  grain  must  be  the  same  as  the  case 
strains  at  that  surface.  So,  combining  Equations  (2a),  (2c),  (4a)  and 
(4c),  we  have 

[cJ.  - v (d  + d )]/E  -€*-€«-  R P /(E  A ) 

L 0p  P rp  zp  p 0p  0c  e e'  N c o' 


(6a) 
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and 


- V]/ep 


€ 

zp 


€ 

zc 


-V  R P /(E  X ), 
c e e'  c c 


(6b) 


where  the  p subscript  refers  to  the  propellant, 
and  (6b)  for  the  axial  propellant  stress  (tf  ) 
interface,  one  obtains 


Solving  Equations  (6a) 
at  the  propellant- case 


( 1 ■=  [cJa  - v 6 - E R P /(E  X )]/v  (7a) 

zp  6p  p rp  p e e'  ' c c/J/  p ' ' 

and 

a = v d +Vtf0  - v E R P /(E  X ).  (7b) 

zp  P rp  p 0p  c p e e'  ' c c'  w 

Equating  (7a)  and  (7b)  yields 

-v  cf  -E  R P /(E  X )]/v  - v (tf  )-v  E R P /(E  X ) (8) 

L 0p  p rp  p e e'  ' c c'  J/  p p'  rp  6p'  c p e e'  ' c c'  ' 

which  reduces  to 

( 1 - v 2)  - V (1  + V ) = E R P (1  - v v )/(E  X ) . (9) 

0p'  p ' p rp'  p'  pee'  c p"  ' c c 

The  radial  and  tangential  stresses  in  the  propellant  (c?^  and 
respectively)  may  also  be  determined  from  the  Lame  equations  for  a c.p. 
grain2  : 

drp  " ^mSPi  " Pe  + ^(Pe  “ VIA1  ~ ni2)  (10a) 

and 

^ep  - - Pe  - H2(Pe  - P^l/U  - in2),  (10b) 

where  m is  the  ratio  of  the  Inside  radius  to  the  outside  radius,  p is 
the  ratio  of  the  inside  radius  to  any  other  radius  (See  Figure  l),  P^ 
is  the  pressure  at  the  bore  of  the  grain  and  P is  the  radial  stress 
at  the  outside  radius  of  the  grain. 
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Evaluating  Equations  (10a)  and  (10b)  at  the  outside  radius  of  the 
grain  (at  the  case  wall  where  p ■=  R /R  = m),  one  obtains 

X 6 


<J  - [m^,  - P + m2(P  - P,)]/(l  - ms) 
rp  Li  e e i ' 


•=  -P±  (11a) 

and 

- Pe  - m2(Pe  - Pi)]/(1  - m2) 

•=  Czm2?^  - P (l  + m2)]/(l  - m2).  (lib) 

Substituting  Equations  (lla)  and  (ilb)  into  Equation  (9)  yields 

[an**.  - P (1  + m2)](l  -v  a)/(l  - m2)  + V P (l  + V ) - 
x e p P e P 

E R P (1  - V v )/(E  X ),  (12) 

pee  c p c c 

which  may  be  written  in  terms  of  the  radial  stress  (Pfi)  at  the  propel- 
lant-case interface  as 


P /P, 

e'  i 


(2m2)(l  - v 2)/(l  - m2) 


E R (l-v  v )/(E  X )-v  (1+v  )+(l-v  2)(l+m2)/(l-m2) 
p ev  c y'  v c c'  pv  p p /v  ' 


(13) 


In  Equation  (13)»  P /P.  is  referred  to  as  the  pressure  ratio,  where 

P may  be  considered  as  a fictitious  pressure  applied  to  the  outside 

surface  of  the  grain  and  equal  in  magnitude  to  the  radial  stress  at 

that  surface.  With  an  assumed  uniform  internal  pressure,  P will  be  uni- 

© 

form  along  the  axis  of  the  grain. 

To  determine  the  bore  stresses  of  the  propellant,  Equations  (10a) 
and  (10b)  are  evaluated  at  the  bore  of  the  grain  (where  p = l)  to  give 


14 

c’rpi  - Cn^i  - Pe  + (Pe  - Pi)]/(1  - m2) 

" "pi  (14a) 

and 

% = ^2pi  “ pe  - (pe  ' Pi)lA1  - mS) 

"O^i  + m2)  - 2Pe]/(l  - m2).  (14b) 


The  external  pressure  on  the  grain  (P  ) varies  only  slightly  from 
the  internal  pressure  (P  ) and  is  constant  along  the  length  of  the 
grain  for  an  assumed  uniform  internal  pressure.  That  is,  the  pressure 
ratio  defined  by  Equation  (13 ) is  approximately  one  and  is  not  a func- 
tion of  axial  location  in  the  grain.  Neglecting  end  constraints,  the 
end  loading  will  be  approximately  equal  in  magnitude  to  the  internal 
pressure.  In  the  actual  rocket  motor,  even  without  end  constraints, 
there  will  be  shear  stresses  at  the  propellant-case  interface  and  the 
axial  propellant  stress  will  be  a function  of  axial  location.  For  the 
current  analysis,  these  shear  stresses  are  not  considered  and  the  grain 
is  under  an  essentially  uniform  pressure  loading  as  shown  in  Figure  5» 

The  axial  propellant  stress  (tf  ) is  then  approximately  equal  to  the  in- 

zp 

temal  pressure  and  is  assumed  constant  across  the  radius  of  the  grain 
and  along  the  axis  of  the  grain.  For  the  case-bonded  grain,  the  pro- 
pellant and  case  strains  were  matched.  For  consistency,  the  axial  pro- 
pellant stress  is  evaluated  at  the  propellant-case  interface  from  Equa- 
tion (7b)  i 


<J 

zp 


v a - verp  /(e  x ). 

p rp  p Op.  c p e e'  ' c c' 


(15) 


Equation  (15)  is  a function  of  the  geometry,  physical  properties  of  the 
propellant  and  motor  case,  and  the  pressure  ratio.  In  addition,  the  ax- 
ial stress  in  the  propellant  is  independent  of  axial  position.  Substi- 
tuting the  radial  and  tangential  propellant  stresses  at  the  case  wall 
given  by  Equations  (lla)  and  (lib)  into  Equation  (15)  gives 

d - -V  P + V [2m2P  - P (1  + mE)]/(l  - mE)  - V E R P /(E  4 ).  (16) 
zp  p e p i e'  J c p e e c c 

The  bore  strains  of  the  propellant  can  be  determined  from  Hooke's 
Law,  Equations  (2a)  through  (2c),  in  terms  of  the  bore  stresses  and  the 
physical  properties  of  the  propellant  and  are 


' - V (<j  + d ), 

©Pi  p lr>i  ZP 

(2a) 

rpi  pv  Bpj^  zp' 

(2b) 

- v (cl  + a ). 
zp  pv  epj^  rpi 

(2c) 

Substituting  Equations  (14a)  and  (14b)  into  Equation  (2a)  gives 

Vepj  ‘ Lf1(i+»a)-a'e]/(i-'"a)+V1+''p2pe  - 

Vp^Spi"Pe(l+m8)^(l""8)+vpvcEpVe/(VC)‘  (l?j 

After  grouping  coefficients  of  and  Pe,  Equation  (17)  may  be  written 


E *=  P,C(l+m2)/(l-m2)+v  -(Em^  2)/(l-m2)]  + 

P 6Pj_  1 P P 

Pe{-2/(l-m2>VpCl+(l-fm2)/(l-m2)>VpVcEpRe/(EcAc)]. 


(18) 


Finally,  the  bore  strain  in  terms  of  the  pressure  ratio,  geometry  and 
physical  properties  is  given  by 


*=  (P^/Ep){[l+ma^v^(l-m2-2m2Vp2)3/(l-ms)  + 


(Pe/Pi)[2(vp2-i)/(l-m2)+VpVcEpRe/(EcXc)]].  (l 

The  radial  and  axial  bore  strains  are  determined  in  a similar 
fashion.  Substituting  Equations  (l4a),  (14b)  and  (l 6)  into  Equation 
(2b)  yields 

E £ = -P.-V  [P, (l+m2)-2P  ]/(l-m2)+v  ap  - 

p rp^  i p i 7 e J/  v p e 

V 2[2m2P.-P  (i+m2)]/(l-mE)+(v  v E R P )/(E  X ),  (2 

pu  le  ' pope  e ' ' c c 


which  reduces  to 


E € = P<L-1_V  (l+m2)/(l-m2)-2m2v  2/(l-m2)"]  + 

p rp.  i p ' p J 


Pe[2vp/(l-m2)+Vp^Vp2(l+ms)/(l-m2)+VpvcEpRe/(EcXc)].  (2l) 


In  final  form,  the  radial  bore  strain  is 


- (Pi/Ep)(-[vp(l+m2+2m2Vp)/(l-m2)  + l]  + 


(Pe/Pi)C2Vp(l+Vp)/(l-m2)+vpvcEpRe/(EcXc)]3.  (22) 


pep 


Substituting  Equations  (14a),  (14b)  and  (16)  into  Equation  (2c) 


gives 


" -VpPe  + vpt2»2I“1-Pe(l+»2)]/(l-»s)  - 

VcVeV(EcXc)  + Vl  * vIP1(l+»")-»e]/(l-»").  (23) 
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After  grouping  coefficients  of  P.  and  P in  Equation  (23)1  we  have 

1 G 

EpG^p  = P^[2m2Vp/(l-m2)  + v^l-(l+m2)/(l-m2)^]}  + 

Pe{2Vp/(l-m2)-Vifl+(Hm2)/(l-m2)]-VcEpRe/(EcXc)}.  (24) 

Equation  (24)  then  simplifies  to  yield  the  axial  bore  strain  of  the 
propellant : 

e*Pl  • (25) 

Note  that  by  expressing  the  axial  strain  in  the  above  fashion,  it 
is  not  necessary  to  know  the  radial  stress  at  the  propellant-case  inter- 
face (P  ).  Instead,  the  internal  pressure  (combustion  chamber  pressure) 
e 

and  the  pressure  ratio  are  required.  The  pressure  ratio  is  a function 
only  of  the  geometry  and  physical  properties  of  the  propellant  and 
motor  case,  so  the  axial  propellant  strain  is  independent  of  axial  posi- 
tion in  the  grain. 

Equations  (19)>  (22)  and  (25)  give  the  tangential,  radial  and 
axial  strains,  respectively,  at  the  bore  of  the  propellant  grain  sub- 
ject to  the  constraints  and  assumptions  of  this  analysis.  Each  of  these 
three  relations  is  independent  of  axial  position  and  is  a function  only 
of  the  geometry,  physical  properties  of  the  propellant  and  motor  case 
and  the  chamber  pressure  of  the  rocket  motor. 


III.  COMPARISON  WITH  THE  FINITE  ELEMENT  ANALYSIS 

The  results  of  this  simplified,  analysis  of  propellant  strains 
are  summarized  in  Figures  6 through  11.  The  strain  of  major  impor- 
tance to  propellant  deformation  is  the  tangential  strain  as  was  men- 
tioned in  the  Analysis  section  of  this  paper,  however,  the  radial  and 
axial  strain  results  are  also  shown.  It  is  necessary  to  examine  each 
of  the  three  strain  results  so  that  a valid  conclusion  may  be  drawn 
about  the  entire  analysis. 

Each  of  the  six  plots  is  shown  in  non-dimensional  form  with  the 
strain  plotted  as  a function  of  axial  location  along  the  propellant 
grain.  The  solid  lines  indicate  the  theory  developed  in  this  paper  and 

4 

the  symbols  represent  the  finite  element  solution  of  Brisbane  . The 
propellant  grains  used  for  these  results  are  140  inch  diameter  c.p. 
grains  with  an  axial  length  of  J00  inches  and  initial  web  thicknesses 
of  40,  34,  28,  22,  16  and  10  inches.  The  results  shown  are  for  these 
configurations,  but  they  may  be  applied  to  any  motor  with  the  same  non- 
dimensional  parameters.  Each  of  the  three  strains  has  six  plots  asso- 
ciated with  it.  Each  of  these  six  plots  represents  a differer t Initial 
web  thickness  (TAU)  which  is  expressed  in  the  non-dimensional  form  of 
web  thickness  divided  by  the  outside  diameter  of  the  grain  (DO).  The 
independent  variable  for  each  plot  is  the  ratio  of  axial  location  (z) 
to  the  outside  diameter  of  the  grain.  The  axial  location  on  the  grain 
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Figure  10,  Comparison  of  the  Finite  Element  Theory  and  Simplified 
Analysis  Results  for  Axial  Bore  Strain. 
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may  be  measured  from  either  the  head  or  aft  end  of  the  grain  since  the 
propellant  deformation  is  symmetric  about  the  midpoint  of  the  axial 
length  of  the  grain.  Each  of  the  six  figures  shows  the  strains  pre- 
dicted by  this  theory  to  be  independent  of  axial  location,  while  the 
results  of  Brisbane  are  definite  functions  of  axial  position. 

The  tangential  strain  plots  (Figures  6 and  7)  show  results  for 
the  six  configurations  used  in  this  analysis.  The  worst  comparison 
occurs  for  the  largest  web  to  outside  diameter  ratio  of  40/140.  For 
this  example,  the  simplified  strain  analysis  gives  a tangential  strain 
that  is  greater  than,  but  within  10^  of,  the  finite  element  results 

over  approximately  one  quarter  of  the  length  of  the  grain.  The  results 

/ 

improve  dramatically  for  the  smaller  web  to  outside  diameter  ratios. 

At  the  smallest  web  thickness  (10  inches) , the  strain  predicted  by  this 
analysis  is  less  than  O.ffl  lower  than  the  results  from  the  finite  element 
solution  over  65^  of  the  length  of  the  grain.  For  this  case,  the  tan- 
gential strain  at  the  ends  of  the  grain  (where  the  largest  difference 
is  expected  due  to  the  uniform  end  loading  assumed  in  this  analysis)  is 
still  less  than  6%  greater  than  the  "exact"  strain. 

The  dashed  lines  on  Figures  6 and  7 represent  the  results  ob- 

2 

tained  for  the  tangential  propellant  strain  from  Vandenkerckhove ' s 
analysis  for  the  rigid  motor  case.  Vandenkerckhove  expressed  his  re- 
sults in  terms  of  the  radial  and  tangential  propellant  stresses.  His 
analysis  is  extended  in  the  same  manner  as  this  analysis  to  obtain  the 
tangential  bore  strain  of  the  propellant.  The  motor  case  and  propellant 
strains  are  matched  at  the  propellant-case  interface  and  are  zero  for 
the  rigid  motor  case.  A uniform  end  loading  equal  in  magnitude  to  the 


! 


r 


I 

I 

I 


I 

I 

? 

F 

r 


27 

internal  pressure  (P^)  is  assumed  (i.e.  Vandenkerckhove  modeled  the 
grain  as  a long,  uniform  thick  walled  cylinder  with  no  end  loading  or 
end  constraints)  and  the  tangential  strain  is  given  as  a function  of 
the  pressure  ratio  and  the  physical  properties  of  the  propellant.  The 
tangential  strain  predicted  by  the  rigid  motor  case  analysis  is  not , 
however,  a function  of  the  physical  properties  of  the  motor  case  as 
is  true  for  the  flexible  motor  case  analysis. 

As  seen  in  Figures  6 and  7,  the  tangential  strains  calculated  from 
the  rigid  motor  case  analysis  are  generally  less  than  the  finite  element 
results,  while  the  flexible  motor  case  results  are  larger.  The  tangen- 
tial strains  predicted  by  the  current  theory  are  closer  to  the  finite 
element  results  than  are  the  rigid  case  results  for  the  middle  portion 
of  the  grain  in  each  of  the  six  configurations.  For  the  larger  webs 
(Figure  6) , the  rigid  case  results  are  better  than  the  flexible  case 
results  near  the  ends  of  the  grain,  but  this  difference  is  less  pro- 
nounced for  the  smaller  webs.  In  general,  except  for  locations  very 
close  to  the  ends  of  the  grain,  the  flexible  motor  case  predictions  of 
this  analysis  yield  a better  approximation  to  the  finite  element  results. 

The  results  for  radial  and  axial  strains  are  not  as  promising  as 
those  for  the  tangential  strain.  The  axial  propellant  strains  predic- 
ted by  the  present  theory  are  close  to  the  strains  predicted  by  the 
finite  element  solution  for  the  smallest  webs  (Figure  ll),  but  in  gen- 
eral, the  results  of  this  analysis  must  be  deemed  invalid  for  both  radial 
and  axial  strains.  In  particular,  the  radial  strains  (Figures  8 and  9) 
predicted  by  this  simplified  analysis  are  greatly  in  error  when  compared 

« 

to  the  results  obtained  from  the  finite  element  method.  This  does  not  * 
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invalidate  the  analysis  for  application  to  the  problem  at  hand,  since 
the  results  for  tangential  strain,  the  only  parameter  needed  for  the 
ballistic  analysis,  are  sufficiently  close  to  the  "exact"  strains  that 
they  may  be  used  to  include  propellant  deformation  effects  in  motor  per- 
formance predictions. 

The  advantage  of  this  analysis,  at  least  for  tangential  strain,  lies 
in  the  simplicity  of  the  calculations  required.  Whereas  the  finite  ele- 
ment technique  gives  more  accurate  results,  it  requires  lengthy  com- 
puter time  and  storage  to  complete.  The  tangential  strain  computed  from 
the  equations  derived  in  this  analysis  requires  only  a knowledge  of  the 
geometry  and  physical  properties  of  the  motor  and  is  done  in  fractions 
of  a second  on  the  computer. 

In  summary,  this  simplified  analysis  for  propellant  strains  yields 
good  results  for  tangential  strain  calculations  and  is  directly  applica- 
ble as  a first  approach  to  including  deformation  effects  in  rocket  mo4  or 
performance  predictions.  Since  the  propellant  deformation  effects  re- 
quire a tangential  bore  strain  calculation  after  each  propellant  incre- 
ment is  burned,  the  disadvantage  of  the  finite  element  technique  is  obvi- 
ous. The  entire  analysis  would  have  to  be  repeated  at  each  bore  diameter 
and  each  such  analysis  requires  substantial  time  to  complcle.  TTie  tan- 
gential strain  calculations  from  this  simplified  analysis  also  have  to 
be  repeated  at  each  bore  diameter,  but  they  require  only  the  evaluation 
of  one  relatively  simple  equation.  This  analysis  thus  provides  a fast, 
simple  technique  for  calculating  tangential  strain  which  may  be  easily 
included  in  existing  software  computer  programs  in  order  to  evaluate 
the  effect  of  propellant  deformations  in  rocket  motors. 


IV.  APPLICATIONS  TO  THE  INTERNAL  BALLISTICS  OF  .SOLID  ROCKET  MOTORS 


I 
I 
I 
I 

| As  was  discussed  in  the  Introduction  section  of  this  thesis,  the 

- propellant  deformations  in  solid  rocket  motors  can  greatly  affect  motor 

performance  and  performance  predictions.  A detailed  analysis  of  the  ef- 
I feet  of  grain  deformation  on  the  internal  ballistics  of  solid  rocket 

motors  is  presented  in  Reference  5« 

I The  basic  hypothesis  upon  which  the  grain  deformation  effects  on 

. internal  ballistics  is  based  is  that,  at  constant  pressure,  the  regres- 

1 sion  rate  of  the  propellant  burning  surface  is  independent  of  the  strain 

I in  the  soil'd  propellant  underneath  the  burning  layer.  In  addition,  when 

Poisson's  ratio  is  approximately  0.5  (as  is  true  for  most  composite 
I solid  propellants)  it  can  be  shown  that  the  density  of  the  unheated  pro- 

■ pellant  is  essentially  independent  of  strain5 . The  unheated  propellant 

is  that  portion  of  the  grain  which  lies  under  the  burning  surface  of  the 
1 grain.  This  thin  layer  of  burning  propellant  is  of  unknown  composition 

(liquid  and/or  solid)  and  the  physical  properties  of  the  propellant  in 

I 5 

I this  zone  are  degraded  to  a state  where  they  cannot  support  shear  stress 

, Thus,  the  density  and  burning  rate  of  the  propellant  within  the  burning 

perimeter  are  considered  independent  of  the  strain  in  the  propellant  be- 
] neath  the  combustion  zone. 

The  internal  ballistics  of  the  SRM  are  modified  for  propellant  def- 

i 
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ormation  effects  by  two  changes.  The  mass  generated  equation  and  the 

time  Interval  required  to  burn  a given  increment  of  propellant  are  both 
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modified.  The  following  discussion  is  taken  from  the  work  of  Sforzini 
and  Footer  and  is  only  summarised.  The  reader  is  referred  to  k. Terence 
6 for  a detailed  discussion  of  the  SRM  performance  prediction  program. 
Only  the  analysis  of  the  deformation  effects  on  internal  ballistics  are 
considered  here.  The  program  modifications  necessary  to  include  these 
effects  are  explained  in  Reference  5* 

The  mass  generated  pr,r  unit  time  per  unit  length  at  the  burning  sur- 
face of  the  propellant  (the  control  surface  is  chosen  to  lie  just  be- 
neath the  layer  of  degraded  propellant)  would  normally  be  written  as 


P r 
o po  o 


(26) 


where  is  the  burning  perimeter  of  the  grain,  is  the  density  of 
the  propellant  and  r is  the  burning  rate  of  the  propellant.  The  sub- 
script o refers  to  the  undeformed  and  unheated  propellant.  The  defor- 
mation modification  requires  that  Equation  (26)  be  written  as 

mG/L  - S0(l  + ee)pporc/[[l  - P(1  - 2v)/Ell  + a(TGR  - TREF)]f  ( 2 7) 

where  C is  the  tangential  strain  in  the  propellant  just  beneath  the  burn- 
ing surface  (calculated  from  the  tangential  strain  relation  derived  in 
the  Analysis  section  of  this  paper),  P is  the  chamber  pressure,  E is  the 
elastic  modulus  of  the  propellant  and  v is  Poisson's  ratio  for  the  pro- 
pellant. The  term  [l+aO^p-T^,)]  is  a further  modification  for  the 
thermal  expansion  effects  in  the  grain.  The  linear  coefficient  of  therm- 
al expansion  of  the  propellant  is  a,  and  TRgp  and  TGR  are  reference  and 
actual  temperatures,  respectively.  Equation  ( 2 ?)  is  presented  exactly 
as  found  in  Reference  5*  The  thermal  effect  is  a very  mall  one  and  is 
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not  directly  connected  with  the  propellant  deformation  due  to  pres- 
surization. Also,  it  is  assumed  that  thermal  and  pressure  loadings 
have  negligible  effects  on  changes  in  the  grain  length  (See  Reference  5) • 
The  time  interval  (At)  required  to  bum  a propellant  increment  (Ay) 
normal  to  the  burning  perimeter  is  calculated  from 

At  = Ay/rQ.  (28) 

The  effect  of  propellant  deformation  due  to  pressurization  and  heating 
will  modify  the  burning  rate  as 

ro  = rc(l  + eeV{[l  - Hi  - 2v)/eIi  + cc(tgr  - tref)]}3  (29) 

where  r is  the  propellant  burning  rate  determined  from  ballistic  test 
c 

motors  or  strand  burners.  The  time  interval  of  Equation  (28)  is  thus 
modified  to  include  the  effect  of  the  deformed  propellant  by  using 
Equation  (29)  for  the  burning  rate. 

Using  the  SRM  performance  prediction,  program  of  Reference  5i  the 
effect  of  the  propellant  deformations  due  to  pressure  loading  are  ex- 
amined for  two  motors.  Detailed  descriptions  of  the  motors  and  the  in- 
put parameters  for  the  deformation  analysis  are  given  in  Reference  5« 

The  thermal  expansion  effects  are  not  present  for  either  of  these  test 
cases  since  the  reference  and  actual  temperatures  are  the  same.  The 
performance  program  is  used  with  and  without  the  propellant  deformation 
modifications  and  the  results  are  compared  with  experimental  data  on  the 
two  motors.  These  results  are  shown  in  Figures  12  and  13. 
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Figure  12  shows  a plot  of  vacuum  thrust  versus  time  for  a Castor 
motor  developed  by  the  Thiokol  Chemical  Corporation  (motor  designation 
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TX354-5)*  This  motor  has  a totally  c.p.  slotted  grain  and  is  thus 
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especially  suited  for  a test  case  on  the  results  of  this  analysis.  Fig- 

7 

ure  12  includes  the  experimental  results  obtained  from  Thiokol  and  the 
results  predicted  by  three  predictions  of  the  SRM  performance  prediction 
program.  Two  of  these  predictions  include  the  propellant  deformation 
effects  while  the  third  does  not.  The  two  predictions  that  include  the 
deformation  effects  differ  only  in  the  elastic  modulus  (e)  of  the  pro- 
pellant. The  nominal  modulus  is  813  psi. , but  the  second  computer  run 
was  made  with  a modulus  of  406.5  psi.  to  demonstrate  the  effects  of  the 
propellant  modulus  on  the  performance  of  the  rocket  motor.  These  effects 
may  be  important  for  a particular  application  since  the  elastic  modulus 
may  vary  between  propellant  batches  and  the  actual  modulus  may  be  uncer- 
tain. 

As  can  be  seen  from  Figure  12,  the  two  performance  predictions  that 
include  the  deformation  effects  more  closely  predict  the  actual  motor 
performance.  The  prediction  with  the  low  modulus  serves  to  point  out 
the  definite  effects  of  propellant  modulus  on  the  performance  predictions. 
This  is  evident  since  the  modulus  of  the  propellant  affects  only  the  def- 
ormation modifications  in  the  performance  prediction  program.  As  the 
elastic  modulus  decreases,  the  propellant  will  deform  more  under  pres- 
surization. This  will  result  in  a larger  actual  burning  surface.  Thus, 
the  apparent  burning  rate  will  Increase  resulting  in  a larger  deviation 
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with  deformation  effects  and  the  low  elastic  modulus  (E=406.5  psi.) 
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comes  closer  to  the  motor  performance  than  does  the  prediction  with  the 
deformation  effects  and  the  nominal  modulus  (E=813  psi.)*  In  the  case 
of  two  or  more  SRM's  firing  in  parallel,  variations  in  performance  will 
contribute  to  thrust  imbalance.  The  effect  of  the  propellant  modulus 
on  the  deformation  of  the  grain,  and  thus  on  motor  performance,  can  be 
helpful  when  investigating  such  phenomena  since  the  physical  properties 
of  the  propellant  may  vary  somewhat  between  motors. 

Figure  13  shows  the  predictions  for  the  Titan  III  C/D  SRM  where 
head  end  stagnation  pressure  is  plotted  versus  time.  Again,  no  scale 
factor  was  used  on  the  burning  rate  and  the  performance  prediction  with 
the  deformation  effect  is  closer  to  the  experimental  data  than  the  results 
without  the  deformation  effects.  This  motor  has  a star  grain  segment  at 
the  head,  end  of  the  grain  but  this  segment  is  treated  with  a burning  rate 
that  is  not  modified  due  to  the  propellant  deformation.  The  present 
strain  analysis  is  not  applicable  to  star  grains,  and  for  that  reason, 
no  deformation  effects  are  included  for  this  portion  of  the  grain.  This 
probably  has  a minor  effect  on  the  performance  predictions  for  this 
motor  since  the  star  segment  is  a small  portion  of  the  total  grain  (i.e. 
initial  star  grain  length  of  95  inches  and  an  initial  c.p.  length  of  613 
inches).  For  motors  with  a star  segment  that  comprises  a larger  portion 
of  the  total  grain,  however,  the  performance  predictions  with  propellant 
deformations  would  be  more  greatly  affected.  The  possibility  of  includ- 
ing star  segments  in  future  work  in  this  area  is  discussed  briefly  in  the 
Conclusion  section  of  this  thesis. 
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V.  CONCLUSIONS 

I 

It  has  been  shown  that  the  strain  analysis  of  this  thesis  has 

applications  in  the  performance  predictions  of  rocket  motors.  Addi-  J 

1 

tional  comparisons  are  needed  to  substantiate  the  accuracy  and  use-  1 

fulness  of  this  technique,  however,  it  is  clear  that  the  theory  can  be  j 

helpful  in  future  SRM  design  and  modification  work.  Such  comparisons  1 

should  include  several  different  motors  with  different  grain  config- 
urations. These  grains  should  include  variations  on  the  straight  c.p.  j 

grain  such  as  slots  or  tapers  as  well  as  entirely  different  geometries 
such  as  star  and  wagon  wheel  grains.  In  the  case  of  the  star  or  wagon  ( 

wheel  grains,  this  analysis  will  have  to  be  modified  to  predict  the 
propellant  deformation.  One  possibility  would  be  to  compute  an  ef- 
fective bore  perimeter  for  these  more  complicated  grains,  or  perhaps  an 
entirely  different  approach  will  be  necessary. 

More  comparisons  on  the  straight  c.p.  grains  are  desirable  in 
order  to  verify  that  the  performance  prediction  results  shown  for  the 
two  motors  examined  in  this  paper  are  typical  for  the  general  class  of 
c.p.  grain  motors. 

i 

Future  work  in  this  area  of  propellant  deformation  should  also  ex-  | 

amine  the  effect  of  changing  the  constraints  and  assumptions  used  in 
this  analysis.  For  example,  the  uniform  end  loading  assumed  on  the 
grain  might  better  be  approximated  by  the  average  of  the  internal  pres- 
sure (P^)  and  the  axial  propellant  stress  at  the  propellant-case  inter- 
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face.  Another  possibility  would  be  to  use  a non-zero  value  for  the 
axial  motor  case  stress.  Some  nominal  value  between  zero  and  one  half 
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of  the  hoop  stress  may  be  chosen,  or  a value  may  be  assumed  based  on  the 
characteristics  of  a particular  motor  (i.e.  thrust  to  weight  ratio  or 
case  properties). 

It  is  important  to  remember  that  one  of  the  prime  features  of  this 
analysis  rests  in  the  simplicity  of  its  application  in  performance  pre- 
diction programs.  That  is,  the  tangential  strain  relations  are  inde- 
pendent of  axial  position  and  are  easily  applied  at  any  bore  diameter. 

An  increase  in  accuracy  may  be  obtained  be  using  a non-uniform  end 
loading  or  a position  dependent  axial  case  stress,  but  this  would  en- 
tail more  complex  calculations  to  obtain  the  tangential  bore  strain. 

In  addition  to  improving  motor  performance  predictions,  the  pro- 
pellant deformation  effects  may  also  be  useful  in  the  examination  of 
other  SHM  ballistic  problems.  The  use  of  burning  rate  scale  factors 
may  actually  disguise  some  ballistic  effects.  For  example,  from  Figure 
13,  one  can  see  that  the  areas  under  the  performance  prediction  curves 
and  the  experimental  data  curve  are  not  the  same.  This  suggests  a 
nozzle  erosion  rate  discrepancy  that  would  be  hidden  if  the  performance 
traces  were  matched  by  the  use  of  a scale  factor5 . A refinement  of  the 
current  propellant  deformation  theory  and/or  the  strain  analysis  may 
lead  the  way  in  analyzing  and  correcting  other  performance  prediction 
errors. 
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